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Abstract 

We use a recent result by Cabezas et al. [T] to build up an approx- 
imate solution to the gravitational field created by a rigidly rotating 
polytrope. We solve the linearized Einstein equations inside and out- 
side the surface of zero pressure including second-order corrections due 
to rotational motion to get an asymptotically flat metric in a global 
harmonic coordinate system. We prove that if the metric and their 
first derivatives are continuous on the matching surface up to this or- 
der of approximation, the multipole moments of this metric cannot be 
fitted to those of the Kerr metric. 

1 Introduction 

In this article we explore the possibilities arising from the method introduced 
by Cabezas et al. [1] in a recent paper (hereinafter CMMR) for computing 
approximate solutions of Einstein's equations which can describe the sta- 
tionary axisymmetric gravitational field of a rigidly rotating perfect fluid. 
In their paper, Cabezas et al. study a fluid with a very simple equation 
of state; the mass-energy density is constant. However, they suggest that 
their approximation scheme may also be implemented with more realistic- 
and more complex-equations of state. Here we use their approach with some 
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minor changes to deal with a polytropic fluid. However, we do not carry 
the approximation as far as Cabezas et al. did. The main reason for this is 
that, in our problem, the metric cannot be written using elementary func- 
tions; it involves a few functions defined by non-trivial differential equations. 
Therefore, we stop at the level of the linearized Einstein theory but include 
quadratic effects in the rotational motion of the fluid. We hope, however, to 
work out true non-linear terms for this metric in the future. 

We aim, first, to show how a polytropic fluid can be incorporated into 
the CMMR framework and set up the mechanism to run. Secondly, we use 
the results the approach provides (which do not differ much from those of 
Newtonian theory) to question the Kerr metric as a "good" metric to describe 
the gravitational field outside a rigidly rotating polytrope. 

The existence of a suitable source for the Kerr metric that is a reasonable 
perfect fluid has been discussed at some length. The first result is due to 
Roos [2], who points out that a perfect fluid cannot be ruled out as an interior 
of the Kerr metric by local arguments based on the constraints that matching 
imposes on the interior metric: there always exists an interior solution in a 
neighborhood of the surface of the fluid. However, this argument is insuffi- 
cient for a problem which involves elliptic differential equations. Moreover, 
the proof that Mars and Senovilla [3J and Vera pi] give for the uniqueness 
of the exterior metric surrounding a stationary axisymmetric perfect fluid 
requires not only the Darmois matching conditions but also an asymptoti- 
cally flat exterior metric. All of this strengthens the case for the solution 
to this problem being global, since it implies that there are conditions on 
the metric in domains of the spacetime manifold which are not close to one 
another: regularity conditions on the symmetry axis, matching on the fluid 
surface, and asymptotic flatness at infinity. Even though we do not actually 
obtain an exact solution to all these problems, we deal with all of them and 
solve them in a coherent way since they are all worked out up to the order of 
approximation we are considering. In this sense, the CMMR approach that 
we follow in this paper casts a reasonable doubt on the existence of a normal 
source for the Kerr metric. 

There is a rather heuristic argument due to Wolf and Neugebauer [5] 
which also suggests that the Kerr metric is not suitable to describe the ex- 
terior gravitational field. Although it does not follow the scheme we have 
sketched above, it too questions perfect fluids as sources of the Kerr metric. 
For this reason we feel it is worth mentioning. 

As we have already mentioned, one of the most important aspects of the 
CMMR approach is its global character. This arises from several assump- 
tions, one of the most important of which concerns the coordinates. They 
are harmonics and cover the entire spacetime manifold in such a way that 
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the metric is of class C 1 on the matching surface (Lichnerowicz matching 
conditions [U]). Another important assumption concerns the metric itself. 
It is assumed that the metric can be expanded in a double power series of 
two dimensionless parameters, A and Q, the first of which takes into ac- 
count the weakness of the field, and the second is related to the rotation 
of the fluid (which is taken to be in rigid motion). This general framework 
is implemented by making other complementary assumptions concerning the 
dependence of the metric on the spherical coordinates, the form of the match- 
ing surface, and the expansion of the metric in terms of the parameters A 
and fl We aim to provide just enough information about all these assump- 
tions to allow the reader to follow our paper easily. However, for a more 
detailed explanation we direct the reader to the original CMMR article and 
to a previous paper by Cabezas and Ruiz [7j. 

In the first three sections of the paper we introduce the polytropic prob- 
lem guided by the CMMR approach. Our aim is to set up the notation 
and almost everything that is necessary to understand the approximation 
scheme. We devote Section 4 to solving the interior problem and Section 5 
to matching it to the CMMR exterior solution. The relationship between the 
global approximate metric we build up and the Kerr metric is discussed in 
Section 6. We prove a theorem that excludes the Kerr metric from the set 
of admissible exterior metrics. In the last section we make some comments 
and develop an argument to extend our result to other barotropic equations 
of state. 

2 Density and pressure 

We consider a Papapetrou-type stationary and axisymmetric metric. We 
can therefore choose coordinates {t, r, 9, <p} adapted to the two commuting 
Killing fields defining the symmetry, and to the two-dimensional surfaces 
orthogonal to their orbits. We denote by £ = d t , the time-like Killing field, 
and by £ = d v , the azimuthal space-like field, so the metric can be written 
as follows: 

+ 7 rr u r ® uj r + -ire{u r ® uj 9 + uj% uj r ) + -f 0e uj 9 ® uj 9 , (1) 

where the 7's are functions of r and 9 alone, and uj 1 = dt, uj r = dr, uj 9 = r d6, 
uj 1 ? = r sin 9 dip. 

We want to link this metric to a rigidly rotating perfect fluid. To do so 
we consider an energy-momentum tensor, 

T=(ii + p)u®u + pg, (2) 
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invariant under the two Killing fields, so the functions fi and p (the density 
and pressure of the fluid) are also functions of r and 9 alone, and u, the 
velocity of the fluid, is a linear combination of the two Killing fields, 



where u is a constant and if), 



{^ tt + 2u j tv r sin 9 + u) 2 7 W r 2 sin 2 9 

-1. 



(3) 



(4) 



is a normalization factor such that, g(u,u) - 

In the framework defined by equations ((H), (j2J) and (J3j), the energy- 
momentum conservation law reduces to a couple of first-order linear dif- 
ferential equations [8], 



<9aP - (a* + £>)<9 a ln?/> = (a, 6, . . . = r ,9) . 



(5) 



These can be integrated if the density and pressure of the fluid are related 
by a barotropic equation of state, p = p(fi). For a polytropic fluid, 



p = kjj, 1+ n (n > 0) 



(6) 



and a simple calculation leads to the following expressions for /i and p as 
functions of the potential if>: 



n+l 



- 1 



n+l 



(7) 



We have chosen the integration constant if)^ in such a way that the pressure 
and density of the fluid vanish on the surface 

£:^M) = Ve, ( 8 ) 

thus defining the boundary of the fluid. 

Equations (J71) and (jHJ) play an important part in our approach. Both 
results are exact for a rigidly rotating perfect fluid. 
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3 Weak field approximation 



We look for a solution to the problem set up in the preceding section. We 
assume it takes the form of a metric g(X, fi) depending on two dimensionless 
parameters A and Q having the following properties: 

1. g(X, O) tends to the Minkowski metric if A goes to 0, that is, 

0(0, Q) = r] = -(J® (J + uj r ® u: r + uj% u e + ^® ^ ; (9) 

we use this limit to identify coordinates {t, r, 9, ip} as the standard polar 
coordinates of flat spacetime; 

2. g(X, Q) tends to a static spherically symmetric metric when Q goes to 
0, which is a solution of Einstein's equations for a fluid with the same 
equation of state, say a polytrope, as that of the stationary axisym- 
metric solution. 

The parameter Q must be related to the rotational motion of the fluid. 
Following CMMR, we assume u = A^firjj" 1 , where r is a constant with the 
dimensions of length that may be identified as the radius of the fluid ball 
in the static limit. The other parameter A obviously accounts for the grav- 
itational field strength. Though we cannot yet identify it with any definite 
quantity (but we suppose that it may be proportional to the quotient of the 
"mass" of the fluid and its "radius" ) we are sure about the role it plays in the 
approximation scheme: we expect the metric to behave as follows for small 
values of A: 

7« « -1 + \f tt , 7 iv « \ V2 ttf tip , 7 W w 1 + A/ w , 

7 rr wl + A/ rr , 7re~A/ re , 7ee~l + A/ ee . (10) 

These expressions agree with the two conditions we impose on the metric 
and they also give a more precise meaning to the kind of approximation we 
are proposing. 

The above expansion for the metric in A leads to a similar expansion of the 
energy- momentum tensor. Taking into account ( fTUi) . it can easily be checked 
that all the quantities entering in T except /x and p have non-zero values at 
A = 0: <7 ~ ?7, V> ~ 1 an d u pa —uj 1 (we use the same symbol to denote the 
vector field and the 1-form). However, a coherent perturbation scheme based 
on the parameter A needs an energy-momentum tensor which tends to zero 
with the first power of A. This can only be achieved if the density and the 
pressure have a linear term in A. This does not seem evident. 
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The normalization factor up to first order in A reads: 



^^i + h(f tt + n 2 v 2 sm 2 e) , (ii) 



where r\ = t/tq. Since the constant ipY, is equal to the value of ip on the zero 
pressure surface, we can assume a similar expansion for it in A and Q, so we 
write 



~ 1 + A 



M + n 2 Q - Mqh^J 



'12) 



This may seem rather bizarre, but it is just a way to introduce two con- 
stants M and k (the reason for this choice will be given later). Substituting 
expressions (ITT]) and (|T2"]) into equation ((7|) and expanding in A, we get 



A 



~ (/« - 2M ) + fi 2 Qr/ 2 sin 2 - ~ + M « 



n 



k n (n + l) 7 

= ^<f, (13) 

\n+l 

This shows that we can make \i a quantity of first order in A if we assume 
k oc A 1_ ~; a suitable choice that takes account of physical density units and 
other technical facts related to the matching is 

1 (™+D\l-n 



k n 4irr^ 



A , (15) 



where r s is a new length constant. This equation can actually be a definition 
of the parameter A whenever we can establish a relationship between the 
length scales, r and r s . We shall do this in Sectional 

We now come to the field equations. We do not aim to work out an exact 
solution but rather an approximate one. Following CMMR once more, we 
use the so-called post-Minkowskian approximation scheme (see for instance 
[9] for a general description of the approach). It requires the introduction 
of new coordinates {t, x = r sin # cos y = r sin 6 sin ip, z = cos 9}; stan- 
dard Cartesian coordinates associated with the spherical-type coordinates 
{t, r, 9, (p}, and the quantity 

h a p = g a p - rj afj , (16) 

which is assumed to be of at least first order in A. Here the indexes a, (3 . . . 
stand for the new coordinates, and expression (JT6l) relates the components of 
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the spacetime metric and the Minkowski metric in this system of coordinates. 
That is, (r] a 0) = diag(— 1, 1, 1, 1) and g a p are some combinations of the 7's 
introduced in (TjQ). These assumptions are consistent with those made before 
in terms of the old coordinates (equation ( TTUT) ) . Moreover, we require the 
new Cartesian-like coordinates to be harmonic coordinates. 

Under the assumptions we have been making concerning the metric and 
the coordinates up to now, the linearized Einstein equations imply that h a/3 
must be a solution of the differential system 

Ah 00 = -2±q n , (17) 

' s 

A 3/2 fi 

Ah 0l = 4 — — q n r] sm9 mi, (18) 
Ah l3 = -2^q n 5 tJ , (19) 

'a 

d k (h k „-~hrj kl ,)=0, (20) 

where A = S^didj stands for the standard Laplacian in Cartesian coordi- 
nates, h = rp^hry^, and (ra*) = (— sin<£>, cos<y?, 0). As defined in ffTHj) . q 
depends on f tt , which is not a great problem since hoo ~ Xfu- Interestingly 
the pressure does not contribute to the right-hand side of equations ( TPTj) . 
(I18p . or (fl9l) . It can easily be checked from definition ffl5|) and formula (Tfj 
that p is a second-order quantity in A. 



4 Slow rotation solution 

In this section we work out an approximate metric which describes the ge- 
ometry of the spacetime we are interested in, which involves the rotation 
parameter Q. We do not intend to find an exact solution to the linear prob- 
lem defined by equations ifPTj) to (EOT) , but rather an approximate solution up 
to order Q 2 . The differential equations used to describe this system include 
one that should verify the component h 00 , which seems to be the most dif- 
ficult to solve, since hoo appears on the right-hand side of the equation in a 
non-trivial way. So, first we deal with this equation and later we try to find 
a solution for the whole system. 

Let us set Q = in ( fl3l) . Substituting the result of ( |T7I) and replacing hoo 
by Xfu (so we eliminate any dependence on A) gives the following equation 
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for the time-time component: 



A/« = - 



2 (\ 



rl V2 



v n 

-fu-Moj . 



(21) 



We can assume spherical symmetry in this limit; that is, /tt(r). Then, by 
introducing a new function, 



fu(r) 



( l' I - ) + Mo 



(22) 



and changing the independent variable, r = r s s, we arrive at the Lane-Emden 
equation of Newtonian theory [10J (see also [11] and [12]). 



$" + + $ n = 



(23) 



We need a solution of this equation which is regular at the origin of the 
coordinates, s = 0. We can select it by looking at well-known results of 
classical polytrope theory: the solution of the Lane-Emden equation that 
satisfies the initial data $(0) = 1 and $'(0) = is chosen (see [10J). Hereafter 
we shall identify the symbol $ with that particular solution of the Lane- 
Emden equation. It is important to remember that there are only analytic 
expressions for $ for a few polytropic indexes, n = (constant mass density), 
n — 1 (the Lane-Emden equation is linear), and n = 5 |10| . In all other cases 
$ has to be calculated using numerical integration methods. 

We now have an approximation of ^ioo which is of zeroth order in Q. To 
obtain a second-order correction, we replace h 00 by 



h o{r,6)tt\f tt (r,6) -> 2A [$(s) + M ] + \n 2 f a (r,9) . 



(24) 



Substituting this expression into (fT71) . expanding this equation to include all 
the quadratic terms in Q, and taking into account that $ verifies ( 12H 
get the following linear equation for the new fu" 



we 



n 



2n 



n-l 



-r/ 2 sin 2 9 h M k 

2 3 



(25) 



This does not appear easy to solve. However, since $ does not depend on 
the 9 coordinate, we can look for a solution of the homogeneous equation 
associated with ( 1251) as a function of r times a Legendre polynomial. The 
CMMR approach [1] sets the rule that a term of order Q 2 in the goo component 
of the metric should not depend on a Legendre polynomial of higher than 
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second order. If we accept this rule, the following function may be enough 
for our purposes: 



4r 2 2 
/«M) = ^0o(s) + 7,(W 

r 6 



2M k + 



Ms) + 3^ 



P 2 (cosfl), (26) 



where 4>q(s) and 02 (s) are, respectively, solutions of the following two linear 
differential equations, 



-0' o + n$ 



n-l. 



i" + -0' + n$ 



;n- 1 



0. 



(27) 



Both equations admit a one-parameter family of solutions which are regular 
at s = 0. We refer to these regular solution when we write the symbols 0o 
and 02. 

We now use the approximate solution for /i o given by (1221 . ( BID , and(l26l) 
to expand the right-hand side of (|T8|) . Eliminating terms in f2 3 amounts to 
substituting q by $, and gives us: 



A/i 



0i 



4 



A 3 / 2 Q 



§ n rj sin 6 1 rrij . 



(2f 



According to the CMMR prescription, the solution of this equation that we 
need has two parts: a specific solution of the full equation and a solution of 
the homogeneous equation that depends on the angular coordinates through 
a spherical harmonic vector, that is: 



hoi{r,e,y) = -x 3/2 n 

where j\ is a constant and 



Ji + 4 $(s) 



sin 9 rrii 



Us) = I" T 2 ^(r)dT. 
Jo 



(29) 



(30) 



This solution is regular at s = and it also satisfies the harmonic condition 
(GO]). 

The equations involving the components hij are easier to solve. Compar- 
ing (iT9l with (JT7J) , it is obvious that h 00 5ij is a specific solution of the first 
equation if h 00 is a solution of the second one. Moreover, the couple h 00 and 
so defined always satisfies the harmonic condition (1201 . We may add to 
any solution of the homogeneous equation and the harmonic condition; 
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that is, the terms associated with the constants a , b , a 2 and 6 2 , but not m 2 , 
in equation (33) of CMMR. However, since they are not necessary to match 
this solution to the exterior one, we omit them. 

Finally, bringing together the results obtained in this section and joining 
the components of the approximate interior metric in a tensor expression we 
have: 



9i 



[-1 + 2A$ + 2AM ) T + (1 + 2A$ + 2AM ) T> c 



4r 



+xn 2 



-^0 o + -(l-^)-2M O K 
o 6 

2 + (T 2 + D 2 ) 



(T + D ) 



Ji 



4 $ 



1 

3 




(31) 



where we have introduced the CMMR notations to denote spherical harmonic 
tensors, T = w* ® u>*, D = S^dx*® dx j , T 2 = P 2 (cos9)u) t ® w*, -D 2 = 
P 2 (cos6>) (Jytte*® cfe J ' and Zi = F^cosfl) («* ® + o;^ g> u>*). 



5 Global metric 

The next step is to connect the interior metric (13 lj) . through the zero pres- 
sure surface defined by with an asymptotically flat vacuum metric. To 
accomplish this we use the CMMR exterior metric, 

9 „«(- 1 + 2A^)t„ + (i + 2^)d„ 

+2Afi 2 ^ (T 2 + D 2 ) + 2A 3 / 2 ^4 z i > ( 32 ) 

where we have set the constants appearing in the original CMMR expression 
equal to zero. This does not imply any lack of generality in our approach, 
since including or omitting these terms is just a matter of convenience for 
matching correctly the interior and exterior metrics. 

Since we have an approximate expression for the interior metric, we can 
use it to get a parametric equation for the matching surface (jSj). However, 
it is preferable to do this using the exterior metric, since we assume that the 
metric is continuous over the zero pressure surface, as in CMMR. This leads 
to . 

r 1 + ^^ + Jf Q [m 2 - 1 ^jP 2 (cos 6) , (33) 
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and also to the expression for ife which we have been using ffT2l . We have 
included a new constant, k, which was absent in the CMMR expression 
for the parametric equation of the matching surface. Even though it was 
superfluous there, we need it to solve the matching problem here. 

We understand the matching of interior and exterior metrics in the same 
way as Lichnerowicz does: the metric and its derivatives must be continuous 
on the matching surface [Bj. However, since we do not have an exact solution, 
we require these conditions to be fulfilled up to the same order of approx- 
imation as the metric is a solution of the field equations. This means that 
we evaluate both metrics on the surface defined by (1331) . then we expand the 
result in Q neglecting terms of a higher order than Q 2 , then we develop the 
matching conditions in spherical harmonic tensors, and finally we equate all 
the coefficients of the expansion to zero to get a set of algebraic equations. 
This procedure leads to the following six constraints, 



Q 2 



2r 2 



M + 



+ 



r s 
2r 



+ n 2 

M 2 



r s M 

M 2 + - + —0 2 — 
3 2r. \r s 



2r. 
1 



r 



+ 



0. 




0. 








fro 


r s 







37 \r. 



r 2 M 



- 2MnK 



4 lQri T (To 
3x + 4 Ji - - + — - 

3 rA \r„ 



4$ 



0. 



(34) 



One may expect a large number of constraints, but as is essentially equal 
to h 00 , all the constraints we can derive from the matching of these compo- 
nents are already included in the matching of h 00 . 

There are two types of constant in the system f)34p . Two which are of 
the first type, ro/r s and Mo, come from the static limit; the others, of the 
second type, appear at a level where rotation is taken into account. This 
last class can be taken as pure numbers when we try to solve the matching 
conditions but this is not the case for the first class, which may be linear 
functions of Q 2 (constants are seen as functions of A and Q in the CMMR 
approach). However, in our problem we mean that it is better to consider r s 
as a true constant and to assume that there is no Q 2 term in the expansion 
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of M . This leads to 

• = 0, 1* (35) 

two well-known predictions of Newtonian theory |10j . In order to justify our 
choice: first, recall expression (fl5l) . which becomes an interesting definition 
of the parameter A if r s is not a free constant. The first equation in (131)]) 
ensures this since it permits r /r s to be determined in terms of the zeroes of 
<3?(s). So we can write, r s = ro/so, where sq is a number 0. Another technical 
reason is based on a careful inspection of the first matching constraint in ( 1341) . 
It seems that any dependence on Q we assign to r s may be absorbed into the 
extra constant k. Second, we do not expect any correction to the mass at 
the Newtonian level; this is the meaning of our assumption concerning M , 
the monopole moment of the exterior gravitational field. 

Taking into account ( |35|) . we can solve the matching conditions ( 1341) to 
get approximate expressions for the first multipole moments of the metric, 

M 2 ^\ + h 2 (so), J x ^\m -\i{s ) : (36) 

the constants of the interior metric, 

«» 2 Tir. A" ;d + «l.). (37) 

SqIVIq 6 

and the initial data needed to pick out regular solutions of the differential 
equations ( 1271) . 

s 1 3 5 

0o(so) « j , 2 So0 ' 2(so) + 2 02(So) + 3 ~ ° ' (38) 

Substituting flHB"]) and fISTj) into expressions fl3T|) and (1321 . we get an ap- 
proximate solution of Einstein's equations inside and outside the polytropic 
fluid up to order A 3 / 2 and f2 2 , which is of class C 1 on the surface of the fluid 
up to the same order of approximation if functions 4>o(s) and 02 (s) satisfy 
the conditions 



6 Kerr metric 

The vacuum metric ( |32i) can be used to describe the Kerr metric near in- 
finity by choosing suitable values for the constants M , J%, and M 2 . This 

1 A classical theorem proves that the solution of the Lane-Emden equation correspond- 
ing to the initial data ^(O) = 1 and $(0) = has a first zero in the interval < s < oo if 
< n < 5 PZ3]. 
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is a straightforward consequence of multipole moment theory in harmonic 
coordinates [13J. So, to the extent that a stationary axisymmetric vacuum 
metric can be identified by its multipole moment s J14j . [T5] we can say that 
the metric ( 1321) coincides with the Kerr metric if: □ 

m = Xr M , ma = A 3/2 ftr 2 J x , -ma 2 = XQ 2 r 3 M 2 , (39) 

where m and a are the standard parameters of the Kerr metric in Boyer- 
Lindquist coordinates. However, if we consider the metric ( 1321) to be the 
exterior gravitational field of a polytrope, the matching constraints restrict 
the values that M , Ji and M 2 can take. Therefore, we may ask if we still 
have enough freedom to make the exterior metric into the Kerr metric. 

Solving the first two expressions in ( 1391) for m and a, and substituting the 
results into the third one, we find that the Kerr quadrupole moment, —ma 2 , 
is a quantity of order A 2 not of order A as it appears on right-hand side of 
the equation. This condition cannot be fulfilled unless M 2 is a quantity of 
order A, M 2 « 0. Equations flHUj) and (|38[) transform this condition on M 2 
into two initial data for the function 2 (s) 

02(S O )~-^, 02(SO)~-T^-- (40) 

3 3s 

We must add to them the regularity condition s = 0. This means there are 
too many conditions for a function defined by a linear second-order differen- 
tial equation. Let us prove that it is not possible. 

Theorem. Let $(s) be the solution of the Lane-Emden equation defined 
by the initial data $(0) = 1, $'(0) = 0, and let s be the first zero o/$(s) in 
the interval (0, +oo), that is $(so) — 0. The differential equation 

( s y 2 )' + _ 6 ) 2 = o ; ( 41 ) 

has no smooth solution in the interval (0, so) such that 

2 (O) = 0, 4> 2 { SQ ) = -\, 0' 2 ( So ) = -A (42) 

3 3s 

To say that 2 (s) vanishes at s = is equivalent to saying that 2 (s) 
is a regular solution of equation (f4T]) . This implies 2 (s) ~ c 2 s 2 (c 2 ^ 0) 
in a neighborhood of the origin. Then, to prove the theorem we have to 

2 For a development of Kerr multipole moments that is closer to our point of view see 
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demonstrate that any solution of this kind can take the values (142]) at s - 
Let us first prove the following: 

Lemma. If a solution 2 (s) 7^ of the differential equation (fJ7]j is regular 
at s = 0, then it does not vanish in the interval (0, so). 
Let us consider the differential equation, 

(s^i) ' + (ns 2 ® 71 - 1 - 2) 0x = . (43) 

Its regular solutions at s = can be written in terms of the Lane-Emden 
solution, 0i(s) = ci $'(s) (ci 7^ 0). If c\ > 0, 0i(s) is negative in (0, so) 
because $(s) is a decreasing function of s in that domain. Let us set c\ = 1. 

Multiplying (T4"Tj) by $'(s) and (T4"3"]) by 2 (s), then subtracting them, and 
integrating the result over (0, s), we get 

W{s) = $'(s)0' 2 (s) - &'(s)Ms) = 4 f $ '(r)0 2 (r)rfr . (44) 

s l Jo 

The integral on the right-hand side is negative near s = if c 2 > and 
positive if c 2 < 0. Let us take c 2 > (the argument runs the same for c 2 < 0). 
Since $'(s) is negative, the product under the integral symbol is negative in 
(0, so) unless 2 (s) vanishes. Let us assume 2 (s*) = 0, < s* < so- Then, 
since 0' 2 (s*) < 0, we have W(s*) > 0, but the right-hand side of (jHJ) is still 
negative at sj Therefore, we conclude that all the zeroes of any regular 2 (s) 
should be outside the interval (0, s ). 

Let us return to the proof of the theorem. We introduce a new function, 



y(s) = s G 



'.(s) 2 



'2 



h (s s 



(45) 



If 2 (s) is a solution of ( )4T|) . y(s) is a solution of the first-order differential 
equation 

v 2 

y' + ^ + ns 6 ^- 1 = . (46) 

Clearly y(s) is the same function for all regular solutions of (T4T1) . particularly 
the solution of fpIB"]) defined by the initial data y(0) = 0. The other two 
conditions in ( f42l lead to a further one on y(s): y(so) = 0. The second and 
third terms in equation f T46l) are positive in (0,so), so y(s) is a decreasing 
function there. Also y(s) is negative near s = if y(0) =0, and it is negative 
in (0, So) if it does not run away to minus infinity for an intermediate value 
of s. From the definition of y(s) and the preceding Lemma, we know that 
this does not happen. Therefore y(s) < at s so y(so) 7^ and (42) is not 
fulfilled. 
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7 Comments 



We have given an proof based on a perturbation approach which rejects the 
Kerr metric as the exterior metric of a polytropic fluid. Even though we do 
not go further than the linear approximation, but consider quadratic terms in 
the rotation, the approximation surprisingly seems to be enough to come to a 
conclusion by analyzing the constraints that the Kerr metric imposes on the 
quadrupole moment of the exterior field. To obtain our result we need to deal 
with a global metric written in a global coordinate system, so we understand 
the matching between the interior metric and the exterior metric in the sense 
that the metric components and their first derivatives must be continuous on 
the matching surface. On this basis one can argue against the generality of 
our result. There are still important questions that we cannot answer: Is it 
possible to get a different global metric by using a less restrictive coordinate 
condition?. Would such a metric admit the Kerr metric as a suitable exterior 
metric?. Our global metric does, however, contain the desired number of 
free constants: lu, which takes account of the rotation; and r , which may be 
related to the central density by means of equations i fTBl) , 

MO) ~ ^ k + 2nft 2 0o(O)l , (47) 
47rrQ L J 

since constraint ( 1381) implies that 0o(O) is a function of so- It is widely thought 
that such a two-parameter metric is the general solution to this problem in 
Einstein's theory as well as in Newtonian gravity. 

Let us give an argument which may extend our results to perfect fluids 
with other barotropic equations of state which are not of the polytropic kind. 
The integrability condition of Boyer's equations ([5]) can be fulfilled by setting 
the density \i to be a function of the normalization factor if). Introducing the 
new variable X = ip/ips, we can write fi = f(X), and the pressure is given 
by 

p(X)=X j\- 2 f(x)d X - (48) 

Then X = 1 defines the fluid boundary, p(l) = 0. It is clear that p is 
a function of /i, at least in a local sense, so the fluid admits a barotropic 
equation of state. We have seen (equations ( flTi) and ( |T2i) ) that to expand 
the metric in A implies X m 1 + Xq. Therefore, what is relevant to our 
approach is the form that the function f(X) takes near X = 1. For instance, 
let f(X) = fi s g(X)(X — l) n (n > 0), where g(X) is a well-behaved function 
at X = 1, g(l) = 1, and fi s is a constant we introduce for convenience to 
make g(X) dimensionless. Then we have fi ~ fi s {X — l) n [1 + o(X — 1)], for 
small values of X — 1, and fi ~ fi s X n q n , for the first term in the A expansion 
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of \i. This expression is similar to the one we obtain for the density of a 
polytropic fluid ffT3"j) . though it will lead to a different value of the parameter 
A. Furthermore, we can integrate ( T481) by using the form of f{X) for small 
values of X — 1, and by expanding the result in A to check that the pressure 
is at least one order higher in A than the density. All of this suggests the 
existence of a common metric up to first order in A for all these barotropic 
fluids. 

Finally, let us note that if the density is not zero at X = 1 (n — 0), the 
resultant first-order metric must coincide with the CMMR metric. It was 
shown in the CMMR article that the Kerr metric does not fit the exterior 
metric of a constant density fluid, and neither does it fit the exterior field of 
this other kind of barotropic fluid. 
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